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The equations obeyed by the vacuum expectation value of the Wilson loop of Abelian gauge 
theories are considered from the point of view of the loop-space. An approximative scheme for 
studying these loop-equations for lattice Maxwell theory is presented. The approximation leads to 
a partial difference equation in the area and length variables of the loop, and certain physically 
motivated ansatz is seen to reproduce the mean field results from a geometrical perspective. 



I. INTRODUCTION 
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o 
o 

<N 

■ The loop representation 0,0,0,0,01 is a useful tool for studying non-perturbative features of gauge theories, both 
in the lattice and the continuous frameworks. Loop-space formulations of quantum gravity have also been developed, 
giving rise to a geometrical setting in which some of the long-standing questions concerning the small distance behavior 
of space-time can be properly addressed 0, EH ■ 

Some years a go, the loop representation formulation of Hamiltonian lattice gauge theories was studied for several 
models [Lll0;El0]- The goal of that formulation was to produce loop-dependent Schrodinger equations, and to 
develop loop-based schemes of approximation to study the spectra and the phase-transitions of the theories. The 
purpose of this paper is to draw some preliminary lines for a similar study within the Lagrangian formulation. 
£vq I Concretely, we shall study the loop-equation that obeys the Wilson loop average (W(C)) of Maxwell theory in the 
lattice. We shall be interested in exploring some basic aspects of the loop content of the theory rather than in 
' obtaining accurate predictions for the relevant observables. 

As a warm-up for the lattice study, we first recall briefly some basic facts about loop-equations for Abelian gauge 
theories in the continuum. As it is well known, the loop-equation [TH ITil UtI Hsj (sometimes called Migdal, Polyakov 
or Schwinger-Dyson equation) for (W(C)) can be exactly solved for the free Maxwell theory within the continuous 
'"j"' ■ framework [Igjj. 

Qh! I 11 the non- Abelian cases, loop-equations have been considered mainly within the large N approximation lip . IrH 
Ibsl |2C| . Some recent applications of the loop-equation approach to QCD can be found in references 0, |23, l23j . On 
the other hand, in the lattice, even the Abelian cases are non-trivial (except in two dimensions), due to the presence 
of the phase transition that separates the confining regime (absent in the formulation in the continuum) from the 
weak coupling one. 

The paper is organized as follows. In section two we review the Wilson loop equations for Maxwell and Chern- 
Simons theories in the continuum from a loop-space perspective. In the last section we consider Maxwell theory in 
the lattice. 

II. MAXWELL AND CHERN-SIMONS WILSON-LOOP EQUATIONS IN THE CONTINUUM 

We shall deal with the vacuum expectation value 

(W{C)) ee J DA W(C) expiS, (1) 

of the Wilson loop 

W(C) ee exp(-ie j> dx^A^x)), (2) 
where S is the action functional of an Abelian gauge theory, that we shall first take as the n-dimensional Maxwell one 

S Maxell =~\ f d n X F^. (3) 
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Since SMaxweii is quadratic in the fields, the functional integration in equation can be performed. Instead, we 
are interested in studying the functional differential equation that (W(C)) obeys, i.e., in the Migdal loop-equation of 
the model. To this end, we shall use a few tools of the loop-space formulation of gauge theories 0, 0, 0, IE H3 ■ 

Consider the space of oriented and piecewise continuous curves in R n . We say that curves 7 and 7' are equivalent 
if they share the same form factor T^(x,-f), defined as 

2^,7)= / dy» 6 (n \x - y). (4) 

Every equivalence class defines a path. The composition of curves, together with the equivalence relation stated 
above, defines a group product among paths. It can be shown that this group is Abelian. Now, let us consider 
path-dependent functionals 'I' (7). We introduce the path derivative S^(x), that measures the change in ^(7) when 
an infinitesimal path u is attached to the argument 7 of ^(7) at the point x, up to the first order in the vector 
associated to the path 

*(«.7) = (l + «' i ^(x))*(7). (5) 

We shall also use the loop derivative A^ix) defined as 

*(<7-7) = (l + «7^A^(i))*( 7 ), (6) 

with cr AlI/ being the area enclosed by an infinitesimal loop a attached at the point x. The loop derivative is readily 
seen to be the curl of the path derivative 

A^ v (x) = d ll 5 v {x) - d u 5^{x). (7) 

Also, we have 

6 fi (x)T»(y n ) = 5»5W(x-y), (8) 

as can be readily shown. 

Since (W(C)) is a genuine path-space function [it depends on closed paths C], it makes sense to take its loop 
derivative. Using the last equation, it is easy to see that 



A^(x)(W(C)) = -ie J DAF^ v (x)W(C)exp{iS}. (9) 

Taking the divergence of this expression and recalling the variational principle for the Maxwell equations, we find 

d»A^(x)(W(C)) — e J DA (10) 

where we have also integrated by parts in the functional sense. Using 

S ^^ = - l eg^(x,C)W(C), (11) 

we finally obtain the desired loop-equation 

d»A^(W(C)) = -ie 2 g^T»(x, C)(W{C)). (12) 

It should be stressed that everything in equation (|12() is a true loop-dependent object, in the sense of the equivalence 
classes mentioned at the beginning of this section. Equation (|12f> is a loop-space version of the ordinary differential 
equation 

dF 

a—(x)=xF(x), (13) 
ax 



whose solution is F(x) — exp(|^-). Following this hint, one can readily obtain the solution 



(W(C)) = Kexp 



d n x / d n yT^x,C)D^{x - y)T v {y,C) 



(14) 
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where D tiv {x — y) is the Feynman propagator of the theory. Finally 



(W(C)) = Kexp 



(15) 



as it should be |30| . 

In the preceding discussion, we ignored that, due to gauge invariance, the Wilson loop average given by equation 
(JIJ is ill-defined. To remedy this, one may add the gauge-fixing term — ^(d^A^) 2 to the Maxwell action J5J). In that 
case, equation ()12[1 changes into 



& l ± ltu (x) + -d v d»8 ll {x) 



(W(C)) = -ie 2 9tlv T»(x,C)(W(C)), 



(16) 



but the net result l|15|l does not change. 

To conclude this section, let us briefly consider the Chern-Simons model, a topological theory whose action is 



S[A] = J c^xe^A^A, 



Now the differential equation for the Wilson loop is 

s^ x A vX (x){W(C)) = -ie 2 T»(x,C)(W(C)}. 



(17) 



(18) 



It is interesting to notice that, being the model topological, its associated loop-equation is metric independent. This 
contrasts with what occurs in Maxwell theory. 

As before, the solution to the loop-equation should be an exponential quadratic in the "loop-coordinates" T fJ -(x, C). 
In fact, it is given by 



(W(C)) = exp 



where 



(x - y) x 
I x - y | 3 



(19) 



(20) 



is the Gauss self-linking number of the loop C, which is a well known knot-invariant. Again, this result coincides with 
what is obtained by performing the functional integral Q for the Chern-Simons model |31| . 

This way of calculating (W(C)) can also be used for other Abelian gauge theories, such as the Maxwell-Chern- 
Simons |32| and the Proca model in the Stueckelberg formulation p9| , since the differential loop-equations are similar 
to the Maxwell loop-equation l|12fl . 



III. MAXWELL WILSON-LOOP EQUATION IN THE LATTICE 



Now we turn to consider loop-equations in the lattice. The partition function of compact electrodynamics, which 
is the model we are going to study, is given by 



with 



and 



Z= / (Q^)exp(-S), 
Jo l>0 



p>0 



(21) 



(22) 



(23) 
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The Wilson functional W(C) is defined, as usual, by 

W(C)=Y[exp(i9i) . (24) 
uc 

In these expressions, I, p and C denote links, plaquettes and closed loops respectively, f3 is the inverse of the tempe- 
rature and 81 is the angular variable associated to link I. 

The Schwinger-Dyson equation for the Wilson loop, that corresponds to the lattice version of equation l|12ll . may 
be written as 24] 

13 J2 A ( l >P) < W (P ■ C) > + A (^ °) < W ( C ) >= °> (25) 
v 

where 



A(l,C) = J2h> , (26) 

and 



I'eC 



1 iorl = V 

6 W = { -1 for 1 = 1' (27) 
otherwise. 

By I we mean the opposite of link I, and < > means statistical average. In equation l|25J) the sum runs over all 
the plaquettes of the lattice, regardless of their orientation. The factor A(l,p), however, restricts the sum to the 
plaquettes attached to the link /. The loop product p ■ C must be understood as the yuxtaposition of plaquette p and 
loop C, yielding a new closed loop. This corresponds to the lattice version of the loop-derivative of Gambini-Tri 'as 
dtl given in equation ||BJ|. Unlike non-Abelian loops, U(l) loops commute (i.e., they are not "ordered") and do not 
have a marked starting point. 

Equation (|25|l arises by equating the averages of the Wilson loop for two different values of the phase factor exp iOi 
assigned to a fixed link I of the lattice. This equality is a consequence of the invariance of the measure under U(l) 
transformations. Expression (|25|) . which indeed represents one equation for each link of the lattice, may be summed 
over all the links belonging to the loop C to produce the single equation 

P E *(!>P) < W (P • C ) > < W ( C ) >= °» (28) 
uc p 

where 

Ao = x;E*"' = E* a£ *' ( 29 ) 

UC I'eC i=l 

with Li being the number of links that appear i times in the loop C . For simple loops (i.e., loops without multiple 
links) Ac coincides with the length of C . 

Equation (|28H (or J2HJ0 is a loop-equation, since any reference to the electromagnetic potential has been summed 
when averaging over the different field configurations. Despite its simple appearance (it is linear and first order in 
the lattice "loop derivative") and unlike its continuous counterpart, it is not possible, as far as we know, to solve it 
exactly (except for D = 2). We shall discuss a simple method to extract some information of this equation that relies 
on similar ideas that were applied to the Hamiltonian lattice formulation in the past jTTt IT^ | . 

To characterize completely a loop one needs an infinite amount of variables. An incomplete list of them could 
contain the length, area, number and types of corners, the class of knottiness to which the loop belongs, and quite a 
few other variables. One could select, on physical grounds, some of these variables and ignore the (necessarily infinite 
number of) remaining ones. Another point of view could be to consider, instead, all the loops up to a given size, 
which leads to a linked-cluster approach |12| ■ Here we shall adopt the first point of view, i.e. .a "collective variables" 
approach, which has already been considered to deal with the Z2 Hamiltonian gauge theory [Uli^l- To motivate the 
choice of variables we are going to make, we can reason as follows. 

It is well known that for "big" loops and D > 3, the Wilson functional decreases exponentially with: a) the length, 
for the weak coupling regime (/? = 1/g 2 — > +00) ; b) the area, for the strong coupling one (<? is the coupling constant, 
or the square root of the temperature) [2^, [2(| [23. For D = 2, the theory can be solved exactly, and the area 
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dependence holds for any value of the coupling constant (i.e., the theory presents a single confining phase). For 
D = 3, approximation schemes and numerical simulations yield that the situation seems to be very similar to the 
D > 3 cases, although the exact solution has not been found [IE EE 113 ■ All this suggests that length and area are 
good candidates to describe, in a first approximation, the Wilson loop functional. Moreover, the length variables Li 
already appear explicitly in the loop-equation l|28l) . With this guide in mind, we approximate < W(C) > in the form 

< W{C) >^f(L u L 2 ,- ■■;A 1 ,A 2 , •••) , (30) 

where Ai is the number of plaquettes of multiplicity i appearing in the loop C . 

This guess about the possible relevant loop- variables should be accompanied by some restrictions on the loops that 
we are going to take into account. Since we are not including "corners" in the list of variables, it seems reasonable 
to ask our loops to be composed of lengthy straight paths, which could be formed by multiple links. Furthermore, to 
avoid ambiguities about how the first term in equation 1|28|) modifies the area variables, we must restrict ourselves to 
consider planar loops. 

We are ready to write down equation 1|28|) for f(Li; Ai). The first term in this equation produces several types of 
contributions, depending on how the plaquette p is appended to the loop C. Let us begin with 



+ {2D - S)iLi(f{L x + 3, U - 1, L i+1 + 1; A x + 1) - f{L x + 3, U - 1, L^ x + 1; A x + 1), (31) 

that corresponds to taking into account all the plaquettes that hit the loop C at links of type i, in such a way that these 
plaquettes lie outside the area of the loop. The first term in this expression arises when the links of the plaquette and 
the loop that make contact have the same orientation. In the second term, the links mentioned above have opposite 
orientations. We have omitted the dependence of / in the variables which are not modified by the attachment of the 
plaquette. 

The second contribution is given by 

t£<(/(Li + 3, U - l,L i+1 + l;Ai - l,A i+1 + 1) - /(Li + 3, L. t - 1,L^ + 1; A t - 1, A^ x + 1)), (32) 

and is produced by plaquettes that lie on the loop area. 
Putting all the contributions together, one arrives to 

iLi \ ~f + f{L x + 3, Li - 1, L i+ i + l;Ai- l,A l+1 + 1) - f(L x + 3, L t - 1, + 1; Ai - 1, A t ^ + 1) 

+ {2D - 3)(/(Li + 3, Li - 1, L l+l + 1; A x + 1) - f(L x + 3, L t - 1, L^ x + 1; A x + 1) 1 = 0, (33) 



which is a difference equation in two infinite sets of variables (the A's and the L's). Based on the exponential 
dependence on area and length that the Wilson Loop should present in the strong and weak coupling limits, we 
consider the ansatz 

oo 

f(L i ;A i ) = l[X^Y^, (34) 

i=l 

which, when substituted into (|33|) . produces the more amenable system of algebraic equations 

(2D - 3)XfXr 1 Y 1 (X l+1 - AVx) + XfXf 1 Yr\X i+1 Y i+1 - X^Y^) + % - = 0, i = 1,2, • • • , (35) 

with Xq — Yq — 1. Equation l|35|l can be seen (for each i) as a single equation for two unknown functions: Xi and Yi, 
Yet, it is not simple enough for our purposes. Hence, we try separately with length-independent and area-independent 
solutions. In the first case, we set Xi = 1 Vi, which yields 

p{Y i+l -Yi_i} + iYi=Q . (36) 

Taking into account the condition yo = 1, the solution to this equation is 

YAm ~ aI (2f3) + bK (2P) ' [67) 
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where li and Ki are the modified Bessel functions and a, b are arbitrary constants. Since li(2/3) must be regular in 
the strong-coupling limit (/3 = 0), we must have 6 = 0; hence 

Yi urn- (38) 

On the other hand, the area- independent solution corresponds to setting Yi = IVi. This leads to the equation 



{D-l)pxl{X i¥l -X i - 1 )+iX i = Q . (39) 

As before, the admissible solution is 

Xl ~ I (i{D- 1)0X1) ' 1 j 

that must be understood as follows. For i = 1, eq. I|40|) gives a transcendental equation for X%. Once this is solved, 
eq. H4()(l again allows us to calculate the remaining X's. 

We recognize that the "area solution" l|38|l is the exact solution of the D = 2 theory j2||, whereas the "length 
solution" l|40(l corresponds to the mean-field solution (which is dimension-dependent), with X\ being the "mean link" 
[27L l2Sj. In the first case, one obtains a dimension independent result that, although being the exact solution only 
for D = 2, it is also a better approximation for the Wilson functional than the mean-field solution in the strong 
coupling regime, even when D > 2. In fact, the strong-coupling mean-field solution is X\ — 0, which corresponds 
to a vanishing Wilson loop. Instead, the "area solution" 1|38[) behaves as < W(C) >= {{3) A for /? — ► 0, which is 
precisely what results from strong-coupling expansions in any dimension. On the other hand, it is well known that 
the mean-field approximation provides a qualitatively acceptable description of the weak coupling region for large 
enough dimensions, although it incorrectly predicts phase transitions for D — 2 and D = 3. 

Summarizing, we have sketched an approach for studying lattice loop-equations that is based on previous work 
in the Hamiltonian formulation. This approach, while being relatively simple, could capture some relevant features 
among those that are responsible for the critical behavior of lattice gauge models. A natural question at this point 
is whether or not there is a simple extension of these ideas to the non-Abelian models, where loop-equations are 
considerably more involved. This and other relevant questions are under work. 

This work was supported by DID-USB grant GID-30 and Fonacit grant G-2001000712. 
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